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t^J- _ Abstract 

Nuclear structure functions and cross-sections for heavy flavour production in lepton- 

[~^. ' nucleus collisions are investigated in the low x region accessible now or in the near future. 

The scattering on a heavy nucleus is described by the sum of fan diagrams of BFKL 

pomerons, which is exact in the high-colour limit. The initial condition for the evolution 

at x = 0.01 is taken from a saturation model, which reproduces the experimental data 

_C on the proton. The A dependence of the structure functions is well described by a 

Mh. power factor A a , with a reaching values as low as 1/2 at extremely low x. The total 

cross-sections for heavy flavour production reach values of the order of mb, and the 

corresponding transverse momentum distributions are sizeable up to transverse momenta 



^ 



Ctf 



larger than the initial large scale \ Q 2 + 4m 



1 Introduction 



In the framework of the colour dipole model [1,2] and in the high-colour limit N c — > oo, the 
scattering on a heavy nucleus is exactly described by the sum of fan diagrams constructed 
of BFKL pomerons, each of them splitting into two. Numerical solutions of the resulting 
equation for the colour dipole cross-section on the nucleus [3-9] were first presented in [8], and 
then in [10-14]. The gluon density introduced in [8] revealed a "supersaturation" behaviour, 
tending to zero at any fixed momentum k as the rapidity Y — ► oo. As a function of In A; it 
proved to have a form of a soliton wave moving to the right with a constant velocity as Y 
increases. This latter phenomenon is related to its scaling property: at high enough Y the 
gluon density at fixed impact parameter b proved to be a function of the ratio k/Q s (Y, b) , 
where Q S (Y, b), the position in momentum space of the maximum of the gluon density, can be 
interpreted as a " saturation momentum" growing as a power of energy. Subsequent numerical 
calculations [10,14] confirmed this scaling behaviour, although the gluon density itself is 
defined differently by different authors. 

Obviously partonic densities in the nucleus are not observable themselves but are only 
certain theoretical tools to calculate the observable quantities. Therefore different definitions 

1 On leave of absence from Department of High-Energy Physics, St. Petersburg University, 198504 St. 
Petersburg, Russia. 



of the gluonic density are well admissible, provided they lead to the same observables. As such 
the structure function of the nucleus is the most directly calculable one. It was first calculated 
in [8] from the found solution of the BFKL fan diagram equation for a very wide range of 
energies (up to rapidities ~ 50) , for the nucleus in the form of a spherical well and for a rather 
specific colour distribution inside the nucleon, not adapted to the experimental data. The 
found structure function proved to be rising as Y with rapidity and roughly linearly with Q 2 
in agreement with the theoretical expectations [7,11]. Later calculations with a more realistic 
nuclear density [15] revealed that the nuclear structure function in fact grows as Y 2 , the extra 
growth provided by the contribution of peripheral collisions not damped by the non-linear 
term in the evolution equation. Still the attention in both [8] and [15] was centered at the 
asymptotic behaviour at high energies, insensitive to the choice of the initial conditions and 
governed only by the internal dynamics of the evolution equation. 

However, from the experimental point of view it is desirable to know the structure func- 
tion at values of x accessible now or in the near future (see e.g. [16]), say down to 1CP 7 
corresponding to Y < 16 -f- 17. At such rapidities the influence of the initial condition is 
still rather noticeable so that its choice becomes a matter of importance. Accordingly we 
performed a new run of calculations employing more realistic colour distributions in the pro- 
ton, following from the saturation approach of [17] consistent with DIS data for the proton 
below x = 0.010. In [10] the parton distributions following from this initial condition were 
studied. In this paper we report on the nuclear structure function following from the BFKL 
fan diagram evolution equation with physically supported initial conditions. The solution 
of this equation also allows us to find the gluon density in the nucleus. We use it to study 
another directly measurable quantity: the inclusive probability for heavy flavour production 
off nuclei (by photon-gluon fusion) in photon induced reactions. 

2 The nuclear structure function from the evolution equation 

In this section we collect the formulae necessary to calculate the structure function from the 
solution of the BFKL fan diagram evolution equation. The nuclear structure function Fi can 
be standardly defined via the cross-sections (Tt,l f° r the collision of the transversal (T) or 
longitudinal (L) virtual photon of momentum q, q 2 = —Q 2 , on the nucleus A of momentum 
Ap: 

F 2 (x,Q 2 ) = ^(a T + a L ), (1) 

where x = exp(— Y) = Q 2 j(2q -p). Both cross-sections can be conveniently presented via the 
cross-section adA(Y, r) for the scattering of a colour dipole of the transverse radius r on the 
nucleus: 

°t,l( y ) = J d 2 rp T , L (r)a dA (Y,r), (2) 

where pr,L is a well-known distribution of colour dipoles created by splitting of the incident 
photon into qq pairs: 

p 2 N r 1 r -> 

PT(r) = -£iY. Z } da{[a 2 + (l-a) 2 }e 2 K 2 (er)+m 2 K 2 (er)} (3) 



2 Other approaches [18] which explicitly contain saturation in the form of multiple scattering and also 
describe the experimental data for the proton, could be used, see [19] for a study of saturation in the low 
Q 2 region within this multiple scattering model. We will restrict ourselves to the use of [17] due to technical 
reasons, see Section 3. 



and 

PL(r) = ^Q 2 J2 Z fJ daa 2 (l-a) 2 K 2 (er). (4) 

Here summation goes over flavours, e 2 = Q 2 a(l — a) + m*, and m/ and Zf are respectively 
the mass and electric charge in units of e, of the quark of flavour /. 

The dipole-nucleus cross-section in its turn can be presented as an integral over the impact 
parameter b: 

a dA (Y,r)=2jd 2 b^(Y,r,b), (5) 

where 2$ has a meaning of the cross-section at fixed impact parameter. The evolution 
equation in Y can be most conveniently written for the function 

4>(Y,r,b) = ^<S>(Y,r,b) (6) 

in momentum space. It reads [8] 

d \ 

— +H BFKL j(f)(y,q,b) = -(f) 2 (y,q,b), (7) 

where y = a s N c Y/Tv and Hbfkl is the forward BFKL Hamiltonian. 

Putting (||) into (|5|) and (Q), and passing to the momentum space, we can express the 
cross-sections cjt,l directly via the function cj) or its 2nd derivative in the logarithm of the 
momentum: 



a T,L =-- ~ I d 2 bd 2 q (j)(y,q,b)wT,L(q) = - \ d 2 b—h(y,q,b)p T ,L(q), (8) 



i 2 bd 2 q<f>(y,q,b)w T ,L(q) = - d 2 l 
TT J IT J q* 



where 



and 



p T , L (q) = J d 2 r pr jL (r) (l - e^ r ) , w T , L (q) = J d 2 r r 2 p T , L (r)e^ (9) 



.■2t-2 , ,., d 2 (f>(y,q,b) 



h(y, g, b) = q z V^(y, q, b) = ™' »' ' . (10) 

Up to a trivial factor, function h gives the gluon density in the nucleus [8]: 

dxG(x,k 2 ,b) 2N Cl2 ~ . / 1 , ,\ 

This definition is not unique (see the discussion in [20]) but natural in the sense that the 
structure function can be directly expressed via the gluon density thus defined. We shall see 
later that this is also true for other observables. 

Note that for numerical calculations the first form in (||) is much more convenient, since 
it requires neither numerical differentiation nor taking into account the singularity at q = 0. 
Functions wt,l can be easily found analytically: 



P 2 N r 1 r / \ i 

Mo) = -g^r v * £ z f J da [ (« 2 + (i - «) 2 ) (i 2 / 2 + g2 ) - m }} J ^ e), (12) 



p 2 N r 1 

' V. „9 ^-v „9 / 

T 

where 



J(,,«)- , 2 h^±g + ». (13) 



The longitudinal density wl is given by the same expression without the term with m 2 r and 
with the substitutions 

a 2 + {l-af -»o(l-a), q 2 /2 + e 2 -> -4e 2 . 

Formulae (Pl)-(|l3D allow to find the nuclear structure function provided function cfi is known 
from the evolution equation (0). To solve the latter one should fix the initial conditions for 
the evolution, that is, the (f)(y, q, b) at the starting point of the evolution which we denote as 
V = 0. 

3 Initial conditions for the evolution 

The initial function (j>o(q,b) = 4>(y = 0,q,b) can be expressed via the dipole-nucleus cross- 
section at y = and fixed b: 

/fj2 poo ( j r 

— e^ r Mr, b) = J o -J (gr)$ (r, b), (14) 

where $o{r, b) = <S>(Y = 0, r, b). 

Staying in the BFKL picture, with a fixed small strong coupling constant a s , at Y = 
the scattering of a dipole on the nucleus is described by the Glauber formula with the dipole- 
nucleon cross-section a(r) given by the two-gluon exchange: 

$ (r,b) = l-e-^ AT ^ r \ (15) 

Here 

a(r)=g 4 f d 2 r' G(0,r,r')pJr'), (16) 



p p {r) is the colour dipole density in the proton and G(0, r, r') is the BFKL Green function at 
Y = 0: 



„/ 



G(0, r, r') = ( 1 + In — I , r > ( < \ = max(min){r, r'}. (17) 



T(b) is the nuclear profile function normalized to 1. Obviously the density p p is non- 
perturbative and not known. If, as in [8], to simplify the calculations one chooses p p to 
be a Yukawa distribution 

e -fJ.r 

Pp {r) = a— , (18) 

then one obtains 

-AT(b)a { r) =5 2C-l + 21nf - Ei(-f) (2 + r 2 ) +e _f (l-f)] , r = pr, (19) 

where C is the Euler constant and dimensionless B combines the information about the 
nucleus density, value of a s and normalization constant a in (|l^). At r — ► and 00, (19) 
behaves as r 2 ln(l/r) and lnr respectively. 

However this sort of initial conditions is only true in the rigorous BFKL approach with a 
fixed and very small coupling. In fact, if one calculates the structure function of the proton 
using (|H]), one finds a very strong Q 2 dependence incompatible with the experimental data. 
To be nearer the realistic situation we therefore choose the dipole-nucleon cross-section a(r) 



in (15) to be in agreement with the DIS data at moderately low x, at which one may expect 



the start of the evolution according to Eq. (]?]). We take x = 0.01 as a starting point for the 



evolution and choose cr(r) as parametrized by Golec-Biernat and Wiisthoff [17] to reproduce 
all DIS data on the proton below x = 0.01: 



air 



(r) = a (l - e- f3 ) , r = /3r, (20) 



with do = 29.12 mb and j3 = 0.234 GeV for 4 quark flavours (with masses m u = m^ = m s = 
0.14 GeV and m c = 1.5 GeV, which we will use in all our computations). Nevertheless, these 
values have been obtained from a fit to experimental data with x < 0.01 and we are interested 
only in <r(r) at x = 0.01, where the quoted values produce a cross-section which overestimates 
the experimental data for F 2 of the proton, while the Q 2 evolution at this x is well described. 
Thus, to better agree with the experimental data at x = 0.01 we diminished <jq to a lower 
value 20.80 mb, keeping the value of /?; we have verified that with this choice, Eq. fll5| ) 
reasonably reproduces F 2 in nuclei at this value of x, see e.g. [21] for another use of (|15| ) in 
this context. Apart from being in agreement with the experimental data, the cross-section 
(|20| ) has an advantage of allowing to relatively simply Bessel transform the Glauberized cross- 
section ( |l5| ) to obtain the initial function <j> in the momentum space according to ( |i~4| ) (see 
Appendix). 

4 Heavy flavour production off the nucleus 

The total cross-section for the photoproduction of heavy flavour / off the nucleus is given 
by (Jt at Q 2 = with wt corresponding to the contribution of this flavour in the sum (12). 
The inclusive cross-section for the production of the heavy quark with a given transverse 
momentum I can be obtained from the expression for the heavy quark density (see [10,22]): 

(a 2 + (1 - a) 2 ) e 2 ^Ki(e6i)Ki(e6 2 ) + [4Q 2 a 2 (l - a) 2 + m 2 ] K (eb 1 )K (eb 2 ) 



bib 2 
d 2 k 1 d[xG(x,k,b)] 



]_ _i_ e -ik(bi-b 2 ) _ g-ikbi _ g ikb 



(2vr) 2 A; 2 d 2 bdk 2 ' ' ^ 

Note that, using (6) and (11) and the fact that &(Y, 0, b) = 0, the last line of Eq. (21) can 
be rewritten as 

^ £ o[$(^bi,6) + $(V,b 2 ,&)-$(y,b 1 -b 2 ,6) 

if we neglect the evolution in Y and take into account only the multiple Glauber rescattering 
in the nucleus implied by (15), we recover the result for the quark density in [22] (Eq. (26) 
in this reference). The inclusive cross-section is expressed via the density (21) as 

da f _ 47r 2 a em _^ 2 f ^ b d[xq f (x,l,b)] 



d 2 l ' Q 2 f J 2 ld 2 b 

Evidently it is finite at Q 2 = with only the transverse part surviving in this limit. 

The leptoproduction cross-sections can be found once the cross-sections for the virtual 
photoproduction are known, which correspond to (21) and ( |22| ) at non-zero Q 2 : 



d<j f {e + A^e' + Q + X) a em f 1 ,\ doj^f + A -> Q + X) 

d 2 UQ 2 dz ~7tzQ 2 \ 2 Z ) d 2 l [ ] 

(see [23] for a discussion on the validity of the equivalent photon approximation in this kind 
of computations) . Here z is the fraction of the energy of the incident electron carried by the 



virtual photon: z = q ■ p/(pi ■ p), where Ap is the momentum of the nucleus and pi that of 
the lepton. 

It is remarkable that both total and inclusive cross-sections for the heavy flavour produc- 
tion are directly expressed via the gluon density (O) in our definition. This perfectly fits the 
assumed photon-gluon fusion production mechanism and thus supports our definition of the 
gluon density from a pragmatic point of view. 



Performing part of the integrations and using Eqs. ( |l0[) and (11), we express the quark 
density as 

d[xq f (x,l,b)] _ Q 2 N C f°° dk . 



I da{ [a 2 + (1 - a) 2 ] I T (e, k) + [4Q 2 a 2 (l - a) 2 + m 2 ] I L (e, fc)} , (24) 

where function h is defined by fllOD, 

It = (2e 2 + k 2 ) X £ - e 2 (e 2 + I 2 + k 2 ) X 3 - e 2 f (25) 

and 

I L = (e 2 + l 2 + k 2 ) X 3 +f-2xt, (26) 

with 

^ = ^TT 2 ' X = J(e 2 + (l + k)*) (e 2 + (l- kf) ' (2?) 

These formulae allow to calculate the quark density both at Q 2 = and at finite Q 2 once 
the gluon density in the nucleus (i.e. function h up to a numerical factor factor) is known 
from the solution of the evolution equation (]7|) . 

5 Numerical results 

Starting from the initial function <J)q at x = 0.01 described in Section 3 we solved numerically 
the evolution equation (^) for values of y up to y ma x = Vo + 4.0, where yo is the initial value 
of y. We have taken a s = 0.2, so that yo ~ 0.88 and the minimal value of x is ~ 10 -11 . 
For the nuclear profile function we use the one corresponding to the nuclear density given 
by a 3-parameter Fermi distribution, with the values of the parameters taken from [24]. We 
have taken into account the contribution of charm but neglected that of bottom, which is 
relatively small in the total structure function. The structure function of Pb in this interval 
of x and for different values of Q 2 between 10 and 10 5 (GeV/c) 2 is shown in Fig. |]. One 
observes that it grows roughly as In (1/x) with x, and as Q 2 with Q 2 . Its ^4-dependence can 
be well represented by a power behaviour 

F 2A (x,Q 2 ) = A a ^^F(x,Q 2 ). (28) 

Power a results dependent both on x and Q 2 . It is presented in Fig. ^. One observes that 

at relatively large values of x a is close to unity and depends on Q 2 very weakly. However 

at smaller x the power goes down up to values below 2/3 naively expected from the nuclear 

screening. At a given x it rises with increasing Q 2 , however the general trend of going down 

with 1/x persists at all Q 2 . 

Our results for the total cross-section for the real and virtual photoproduction of charm 

and bottom (with m& = 4.75 GeV) on various nuclei are presented in Figs. |3| and ||. The 

cross-sections are shown as functions of the photon cm. energy W with respect to a nucleon 

in the nucleus. We assumed that 

4m 2 + Q 2 

x = W 2 + Q 2 • (29) 



The minimal value of W corresponds to the maximal value of x = 0.01 for our solution. 
Accordingly it rises with Q 2 , which explains why the curves for different Q 2 start from 
different W. As expected, the cross-sections rise with increasing W and 1/Q 2 ■ Their absolute 
values are rather big, reaching ~ 3 and ~ 1 mb for charm and bottom respectively for Pb 
at W = 1000 GeV and Q 2 = 0. The growth with W in the shown interval W < 10 4 GeV is 
rather close to power-like. However, corresponding to the behaviour of the structure function 
(Fig. [y), at higher W this should change to ~ In W. 

In Figs. H and ^ we show the transverse momentum distribution of the produced charmed 
quark off Pb at Q 2 = and 30 (GeV/c) 2 respectively, and at W = 200, 1000 and 5000 GeV 
(in mb/(GeV/c) 2 ). Here we assumed that 

Ami + Q 2 + 4l 2 

x = W<? • (30) 

Our restriction x < 0.01, combined with discrete values of y at which the gluon density was 
calculated, severely limited the number of points at small W. This explains a somewhat 
irregular form of the curves at W = 200 GeV. One observes that, in accordance with the 
BFKL kinematics, a relatively large part of the charmed quarks are produced with momenta 
of the order or even larger than the initial large scale \/Q 2 + Ami . 

6 Discussion 

We have calculated the nuclear structure functions at low x which follow from the fan dia- 
gram evolution equation in the perturbative QCD, with the initial conditions adjusted to the 
existing experimental data on DIS on the proton at x = 0.01. Knowing the gluon density in 
the nucleus we also calculated the total and inclusive cross-sections for charm and bottom 
production at low x. The results show that the structure functions grow with increasing 1/x 
and Q 2 as In (1/x) and Q 2 respectively. Parametrizing the A dependence as A a we found 
that at very small x, a is going down to values below 1/2. The found transverse distributions 
of heavy quarks exhibit a relatively large contribution from the momenta comparable or even 
larger than the natural scale J Ami + Q 2 . 

Our results have been obtained in the pure BFKL kinematical regime without any addi- 
tional cuts. However they are based on the approximation of large N c . Therefore one should 
not expect that they remain valid at extremely small x when the 1/N 2 corrections growing 
with 1/x may change the found behaviour. 

We can compare our results for the structure functions with the recently obtained in [12]. 
Although the general trend is similar, our structure functions grow with 1/x considerably 
faster. In particular, for Au at Q 2 = 100 (GeV/c) 2 and x = 10~ 7 , and also taking a s = 0.25 
as in [12], we obtain F^a/A ~ 140 as compared to ~ 20 in [12]. The difference may come 
partly from a different initial condition, which in [12] was taken in a simplified form (with a(r) 
roughly proportional to r 2 ) and partly from the fact that in [12] the asymptotical form for 
the photon dipole distributions (||) and (||) at large Qr was used. The latter approximation 
seems of dubious validity to us, since in fact values Qr ~ 1 give the bulk of the contribution 
to the structure function. Besides, in [12] the integration over the dipole dimension r was 
cut at both small and large values, which obviously introduced a dependence on the cut-off 
parameters. 
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7 Appendix. The initial function with the Golec-Biernat- 
Wiisthoff dipole cross-section 

Using Eqs. (|lj|), (|i~5| ) and (|2(]) and passing to the integration over x = j3r we write the initial 
function cfro(q,b) as an integral 

Mq,b) = J o ™Mqx/p)\i-e l J J. (31) 

The dependence on b is contained in the dimensionless factor B: 

B = ^AT(b)a . (32) 

In the following this dependence will be suppressed. We will also denote z = q/ (3. 
We present the exponential as a power series in its argument: 



-B 

1-e 



1 -x 2 

1— e 



E^F (-«-)*. m 



71=1 

so that 

n=l 

where 

In(z)= r ^J ( XZ ) (l - e- x2 Y . (35) 

Jo x V / 

The integral I n can be presented as 



n /"OO ,7™ „ n /-OO Arf , „. 

In(z) = Y,C*(-l) k ^Uxz)e~^ = -Y,C k n {-l) k -JoM 1-e"^ 2 . (36) 
t^n Jo x f~r Jo x v / 



—Uxz)e~ kx2 = -J2r L -'■■ /x llr 

k=o " u x k=0 

The last integral in x can be found analytically. Indeed one has (see [25]) 



o 2^+ 1 a(^+ l '+ 1 )/ 2 r(^ + 1) ^ 2 ' 4a y v ; 

Taking v = one has 

**"«-*(») -^^^^f^i;-^)- (») 



We take the difference of two integrals ( pq ) with different a's, 

f°° dx x» (e~^ 2 - e~ ax2 ) 3 (zx) = 



f^l;i;--l -a-^l\ Fl fall;!;-^ 



ir((/i + l)/2) T^^iF, 

The result is obviously finite at \i = — 1. To find it we put 

// + 1 = 2A -»• 0. 
In this limit, 



7 



~ A ~1-Ahi7, a _A ~l-Alna, r(A) ~ 1/A 



and 



iFi(A;l;t) ~ 1 + A(Ei(t) - C - ln(-t)J , i < 0. 
Putting this into (|39"|) we find, in the limit A — ► 0, 



dx 



-•yx 



J (zx) 



/o £ 

Taking here 7 — ► we find the desired integral: 

roc dx 




Ei 



4a 



|°°^(l-e— 2 )j (^) 




(39) 

(40) 

(41) 
(42) 

(43) 



(44) 



With the help of this formula one can find the integrals I n (z) as a finite sum of exponential 
integrals of different arguments, Eq. (p6|). Putting this into ([34|) gives the initial function fa. 
For realistic values of the parameter B < 3 the convergence of the series in n is very fast, so 
that (J)q can be calculated with very high accuracy at all q. 
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Figure captions 

Fig. 1: The structure function F2 of Pb as a function of x at different Q 2 . 

Fig. 2: The power a of the A-dependence of the nuclear structure functions as a function 

of x at different Q 2 . 

Fig. 3: The total cross-sections for real and virtual photoproduction of charm on Pb (lower 

right plot), Ag (lower left plot), Cu (upper right plot) and Ne (upper left plot). 

Fig. 4: The total cross-sections for real and virtual photoproduction of bottom on Pb (lower 

right plot), Ag (lower left plot), Cu (upper right plot) and Ne (upper left plot). 

Fig. 5: The inclusive cross-section do /d 2 l for the real photoproduction of charmed quarks 

on Pb, as a function of the photon-nucleon cm. energy W. 

Fig. 6: The inclusive cross-section do~/d 2 l for the virtual photoproduction at Q 2 = 30 

(GeV/c) 2 of charmed quarks on Pb, as a function of the photon-nucleon cm. energy W. 
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Figure 2: 
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